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ABSTRACT 


There were two objectives in the experiments performed on the 
barotropic model. One was to correct the lag of the troughs and 
ridges noted in the barotropic model by the Naval Weather Research 
Facility, and the other was to introduce a baroclinic development term. 

The approach to these problems was one which could be incorpo- 
rated into the Fleet Numerical Weather Facility model with ease. 

Included herein are the theoretical development of the numerical 
model, a comparison with the FNWF model and charts illustrative of 
the results with various values of the pertinent parameters for 24 
and 48 hours. 

It was observed that the development term causes selective 
development of lows and highs; but when overdone, spurious short 
waves are formed which, however, do not basically change the major 
wave location. 

It was also observed that the function of latitude used here to 
correct the lag in the horizontal advection of vorticity, resulted in 
greatly improved movement of troughs and ridges from 30 degrees 
to 60 degrees of latitude, but could be further refined to produce even 


better results. 
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no development term. 
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chart V.T. 0000Z 26 Jan 1967, 
term. 


chart V.T. 0000Z 26 Jan 1967, 


no development term with heavy smoother. 
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term, 1/2 of heavy smoother. 


chart V.T,. 0000Z 26 Jan 1967, 
term, 1/4 of heavy smoother. 


1200Z 25 Jan 1967 500-mb analysis 


Page 
Ze 


Zo 


24 


28 
“ae, 


34 


55 


34 


36 


ow) 


39 


40 


41 


43 


44 


lg eae at Page 


16, 48-hour prognostic chart V.T. 1200Z 25 Jan 1967, 45 
no development term with heavy smoother. 


l7. 48-hour prognostic chart V. TigliZ00 Za2bJanwlOom Ab 
1,0 of development term, heavy smoother. 


18. 48-houreprognosticechartpViede W200Z 25 Janeb9Oo#, 48 
no development term 1/2 of heavy smoother. 


19, 48-hour prognostic chart V.T. 1200Z 25 Jan 1967, 49 
1.0 of development term, 1/2 of heavy smoother. 


20. 48-hour prognostic chart V.T. 1200Z 25 Jan 1967, 50 
1.0 of development term, 1/2 of heavy smoother 
twice every 6 hours vice three times. 
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Coriolis parameter 

Relative vorticity 
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Horizontal wind velocity 

Horizontal thermal wind velocity 

Relative thermal vorticity 

Height 

Acceleration of gravity 

Height less standard height 

A scaled parameter 

Map factor 

Distance between grid points 

advection of horizontal vorticity tuning constant 
development term tuning constant 

smoothing factor 
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1. Introduction 

The numerical model for the prognosis of the 500-mb level used 
by the FNWF is based upon the barotropic potential vorticity 
equation [3]. This model has been carefully ''tuned" and engineered 
over a period of several years and is found to produce quite satisfactory 
results. The use of a numerical model derived from a baroclinic 
vorticity equation, however, should improve prognosis by permitting 
development and decay of pressure systems. Due to time limitations, 
it was decided to utilize a quasi-baroclinic model which makes use of 
of the initial (1000-500)-mb thickness anda simplified prognosis 
thereof to compute a Sutcliffe development term in the 500-mb prog- 
nostic equation. No attempt is made to verify the thickness forecasts 
nor the implied 1000-mb prognosis. 

The model developed and tested here is still essentially a one- 
parameter model and constitutes only an attemp to provide for some 
development in the conventional barotropic model, which is the primary 
shortcoming of an otherwise very successful prediction scheme for 
mita-troposphere, 

In addition to the use of a baroclinic term, some experimentation 
was performed on the horizontal vorticity advection term of the 
equation by varying it as a function of latitude. Various degrees of 
smoothing were also applied to the prognosis and the effects on the 


appearance and verification are noted. 
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2, The FNWF model [3] 
The barotropic potential vorticity equation is taken as 
D Ct+T 
cae 
(2.,1) 

where h is the depth of fluid being considered. For application to the 
atmosphere, this layer is considered to be approximately the tropo- 
sphere and to be centered at the 500-mb level. Expanding equation 2. 1 


and multiplying through by h gives 


LED) LDU TLE THEO on 
The last term is assumed to be composed of three components; 
namely, 
a) a mountain term, 
b) a latitude term, 
c) a temperature term. 


Equation 2.2 then becomes 


| ¢ 
aae (48) RR VIVE OS - (HY Vh, 
(1) (2) (3) (*) (=) (2H 3) 
-L(fYVVh- + (ft AVh, =O 
(6) (7) 
Term (2) is the Helmholtz term. It is approximated by 
gil thea CIB =~ poe 
h (H+) 53 rae oO (2. 4) 
where Z is the height of the 500-mb surface, f is a mean value for (f +), 


h is a mean value for h, andyw is a constant determined by experiment- 


ation. 


Term (5) is a mountain term. The approximation 


| 
=H (f+ 2)VVh, 2 at ViINemwy ih (2. 5) 
h 
is made where M is the height of the terrain, and "'a''! is an emperically 
determined constant. With other considerations and simplifications 


this results in 


. = 
(VNA)ZE = —Vi7 (tr Se of M— bEH ) 
ot tT 
(Ze) 
where ''b'' is a tuning constant and H is a (500-1000)-mb thickness. 
The thickness is carried forward during the forecast ag an 
independent field with half of the speed of the 500-mb wind. 
In addition to the above, other engineering and climatological 


modifications are made in the FNWF model but will not be discussed 


here, 
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3. The baroclinic model 

The baroclinic model used here is developed from the horizontal 
advective model [2, or 400) in which temperature changes are assumed 
to result from horizontal advection only. In this model, it is also 
assumed that the vertical variation of the magnitude of the thermal 
wind is the same everywhere and the wind field may be expressed 

Vay + Alp) > (3586 

where W is the vertically integrated mean wind with respect to 
pressure and V.- is the vector difference between the wind at two 
selected isobaric levels. With the geostrophic assumption, Ve is 


the thermal wind. 


By the same assumption 


ye Fs A(t, i) 2) 


Integrating the vorticity equation 


22 ewe ltt) = fF 2 (3.3) 


after substituting from (3.1) and (3. 2) gives 


; ne 
£) FAA ee) Crawnng) (FA %ef)dp 


: (3. 4) 
= $W, 


rY:y-) 
The left side leads to 


2 -V.% (Ff) + PM Vor (3. 5) 
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By the geostrophic approximation 


f = Lv (3. 6B) 


Ve = bkxVH (3. 6C) 
i = Vk (3. 6D) 


substituting (3.6) into (3. 5) 


ret +thevZ- ‘Tol Prf)+ A Lexy. Bev H (3.7) 


leads to 


- wg 2 
22 4 (- 22 29 22,55 22) A) AP 2H 2 ml rf +t Ev H ) (3. 8) 
using 


J(A,8)= G32 — 38-3) 


3.8 hecomes 





2 
7732-4 J(2,0) +A TH L0H) (3.9) 
The right side of (3.4) leads toa Helmholtz term for which the FNWF 


version is utilized. Introducing tuning coefficients K, and K, yields 


to the final operational equation 


(> A)SE = — KT (2,9) “Ky THE 04H) inks 
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where Z is the 500-mb height, H is the (500-1000)-mb thickness, and 
A is the FNWF Helmholtz term. Assuming the thickness changes are 


purely advective gives 


2th, Ws He og SR H2-T(2, H) (3. 11) 


where K, is 0.5, 


V 
is 
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4. The numerical form of the baroclinic equation 

Equations (3.10) and (3.11) are solved over a field composed 
of 63x63 grid points which are a distance d apart. The appropriate 
map factor m for the FNWF 63x63 grid is that for a polar stereo- 
graphic projection (i) . In order to utilize the FNWF subroutine 
library and reduce computation time, fixed point arithmetic was used. 


The constants and corresponding scaling factors are: 


d = 381KM = 3.81 x 10! cm 
m = oe 

GeO essai) cae u 
g = 980 cm sec 


A 17 | 
D=Dx2 ° (D= Z-Zstd) 


A 
ee ee 


~ 9 


? =4 x 2 
A hy, 
= 2 9 % * 
Hy * Sr * (zV'H=t) 
At = 1 hour = 3600secs 


which result in 


(Vem)aban Kan 02eKT(8,5)—Kaxax FPA) (4.1) 
which is the form of equation (3.10) used in the FNWF Helmholtz 
subroutine fr) j 

The thickness is carried forward during the forecast as an 


independent field with half of the speed of the 500-mb wind 


ee 


ot VH = -4-T(2,H) (4. 2) 
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When expressed in finite difference form, (4. 2) becomes 


A A oy 
a he-ax.79Sxam*x J(B,H) (4, 3) 
Both (4.1) and (4.3) are computed on an hourly basis. The first 


time step is produced by the method of forward time differences 
D =D + 4D (4, 4) 


If continued,this method would lead to computational instability. For 
this reason, all other time steps are produced by the method of 
centered time differencing 4 


mt+1l a pe 


D = D (4, 5) 


In the event a smoothing operation is performed ona field, a forward 


time difference is used as the subsequent time step. 
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5. Research procedures 

Comparing equations (2. 5) and (3.10) shows that the FNWF model 
without certain refinements closely approximates the baroclinic model 
when K, ls zero. For this reason, this baroclinic model without the 
development term can be referred to as the quasi-F NWF model; 
however, the FNWF model has additional refinements not included 
here, hence comparisons of verification scores are not strictly valid. 

The data used for the experiments were recent 500-mb and (500- 
1000)-mb thickness fields obtained from FNWF. The actual computing 
was done on the Postgraduate School Computer Facility CDC Model 
1604, 

The results of each experimental run were tested by the pillow 
and RMSE (i) , then analyzed and visually compared with the analysis 
at verification time as well as with other experimental runs. 


Pillow E(A-0y 


x 
| | ee 
RMSE Ss [ (4-1) Pillow Ni 
N=| x 


It is apparent that the RMSE, although a rapid method of 





quantitative comparison of fields, will not show a variation of develop- 
ment of systems. It was therefore necessary to use the visual method. 

The experimental procedure consisted of varying the values of K, 
or K, or both. The values of K, varied from 0 to 2.0 and also 


included variations of the tuning constants with respect to latitude. 
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In addition to the above, the fields were compared after different 


degrees of smoothing were applied. 
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6. Experimental results 

Consider first the advection term only. The speed of movement 
of the troughs and ridges were modified by the use of K, as a function 
of latitude. The need for this modification was made evident by a 
study of the FNWF model conducted by the Naval Weather Research 
Facility, Norfolk, Virginia, which indicated consistently slow move- 
ment of troughs and ridges through certain latitudes, a deficiency that 
appears to characterize barotropic models in general and even some 
baroclinic models. 

After several trials, an emperical correction was made for this 


lag by the employment of a function of latitude (SF) as the factor 


5+ SIN | 
Se Sint 


applied to the horizontal vorticity advection, 

The graph of this function, Fig. 1, shows a maximum of 1.37 
at 21.4 degrees and a value of one at the equator and pole. 

The effect of the use of this sine function can be seen by comparing 
Fig. 2, the 24-hour prognostic chart without the sine function and 
Fig. 3, in which the sine function is applied. 

There is an obvious displacement of all systems while the 
contour shape remains essentially unchanged. 

To test the effectiveness of this function on independent data, a 
test was run on seven consecutive 24-hour periods from 0000Z 24 March 
to 0000Z 27 March. The displacement of troughs and ridges was 


measured in degrees of latitude and an average deviation from the 
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TABLE I 
24-Hr Prog 
AVERAGE DISPLACEMENT FROM 


ANALYSIS 
IN DEGREES OF LATITUDE 


LAT WITH SINE FUNCTION WITHOUT SINE FUNCTION 
TROUGHS RIDGES COMBINED TROUGHS RIDGES COMBINED 


10 -2.72(11) +-2.0(2) -2.62(13) -2.64(11) -2.0(2)  -2. 54(13) 


20 1.43(23) -.36(14)  .76(37) .61(23)  -.93(14) —-. 03(37) 
30 .70(36) -1.91(23)  -.32(59) ~~. 97(36) = -2. 26(23) -1. 47(59) 
40 .51(35) -2.40(30) -.83(65)  -.66(35) -3.57(30) -2.00(65) 
50 -. 35(31)—-- 2. 15(26) -1.16(57) -1.52(31) -3.46(26) -2. 40(57) 
60 -1.66(18)  -.95(20) -1.30(38) -2.38(18) -2.15(20) -2. 26(38) 


(# of samples) 


aS 


ie I oir 
24-Hr Prog 
TEST OF SINE FACTOR 


PILLOW* RMSE* 


Vou NO SF WITH SF FNWF 
00Z 24 Mar a ye: = 94 <1 6104 
12Z 24 Mar =5 2G 25) INGA Bese 11.15. 
00Z 25 Mar aS 90 5. 91 oe Oe We 
12Z 25 Mar 5 Omi 5 watt = Site eS. 
00Z 26 Mar OulS Otel 1:3 = 8 
ZZ 20" Ma = (iy ond LLL Y) =i alls “lace Leb 
0O0Z 27 Mar 18 94 18 2 =25.5 Cm 
*feet 


26 


corresponding analysis was taken for the entire period. 

The results of this test are summarized in Table I. The first 
item observed by this summary was that the movement of troughs 
and ridges was truly in need of an increase at almost all latitudes 
tabulated. It is seen that the use of the sine function produces a 
distinct improvement from 30 to 60 degrees. This is particularly 
evident in the ''combined" column which includes troughs and ridges. 

Table I also shows a possibility of too much correction at 20 
degrees. 

The results at 10 degrees are considered to be inconclusive due 
to the small number of samples and the difficulty of defining troughs 
at this latitude. 

Table II compares the pillow and RMSE for the verification dates 
of the test without the sine function, with the sine function and for 
the FNWF product, The use of the sine function produces little change 
in this statistic. 

Turning now to the development of systems, the problem was to 


define a proper K, term or to determine how much development term 


3 
FO Use. 

The experimental data was selected at random from the FNWF 
files of height fields and corresponding thickness fields. 

The 24-hour prognostic tests were run from 0000Z 25 Jan 1967, 


Fig. 4, to 0000Z 26 Jan 1967, Fig 5. This date was chosen at random 


and, although not the best of all possible situations to illustrate 


development, does show the effects of the several variations. 
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For purposes of comparison, the following factors are considered: 

a) the pillow(P) and RMSE(R) in feet 

b) illustrative charts 

c) the mean of the deviation of selected low and high central 
values in meters from the analysis, together with the standard 
deviation from this mean. 

It will be seen that a low RMSE does not necessarily indicate a 
good verification and therefore must be carefully used in conjunction 
with other evidence before drawing conclusions. 

For information, the FNWF pillow and RMSE for this date 
were 3.4 and 115.9, respectively. The mean of the deviation of lows 
was 7.8 with a range of 92 to -46 and a standard deviation of 47.1, 
and the mean of the deviation of the highs was 11.7 with a range of 
31 to -2 anda standard deviation of 14.0. 

Caution should be used in comparing the FNWF barotropic results 


to the barotropic model used here, where K. is 0, as this model only 


3 
approximates the FNWF model. However, comparison of cases with 
and without the development term should indicate what effect the 
development term would have onthe FNWF model. 

The 500-mb analysis during the forecast period showed little 
significant movement or development with the exception of the low 
pressure system on the northwest coast of the U. S. which moved 


inland and filled, while the low in the Gulf of Alaska deepened by 90 


meters and moved east leaving a trough trailing behind. 


30 
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Case I 
Fig. 6 is the result of using no development term, K, = 0. The 
general position and intensity of pressure systems are good and 


follow the history intensities. The low on the west coast of the U. S, 


fails to fill while the low in the Gulf of Alaska deepens by only 20 


me tense 

mean deviation and range standard deviation from mean 
P -7 lows -5.0(51, -38) 29.6 
R 105 highs 19. 2(28, 9) 6. 8 


K, was next set at 1.0. The prognosis became more unstable 


with each time step and after 13 hours, the FNWF Helmholtz sub 
routine (] would not converge. 

Various other values of K, were tried. It was found that when K, 
is 0.5, the model runs for 18 hours, K, at 0.3 for 23 hours and K, 
equal to or less than 0.2, over 24 hours. 

Since the hourly change of the thickness field is based solely on 
horizontal advection, it is believed that this field becomes less 
representative of the actual field fairly soon and is probably a major 


source of observed instability. For this reason, K, was divided by 


3 
@ every six hours in the following examples unless otherwise 
indicated. Other periodic reductions of K, were tried, however, this 
one gave the best results. 

In addition, the FNWF subroutine SAU fi) is employed on the D 
and H fields every three hours to suppress the generation of short 


Waves. 
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Case Il 

ry eee K, is 0.5. Although the major features remain 
unchanged, compared to Fig. 6, there is some increase in short wave 
development. The mid-latitude lows are from 20 to 40 m deeper and 


highs 10 to 30 m higher. 


mean deviation and range standard deviation 


P -8 lows -8, 3(45, -39) fee Th 
R 115 highs 26. 4(50, 8) 14.9 
Case III 


in Figs, K, is increased to 1.0. Further short-wave development 
is noticed in addition to stronger development of both highs and lows. 

The development is selective, as can be seen by the low north of 
Siberia remaining unchanged while the low in Idaho deepens by 50 m 
and the one in the Gulf of Alaska by 80 m. 

mean deviation and range standard deviation 

Pp -8 lows -16. 8 (2, -43) 152% 
R 128 highs 34. 6(62, 8) (fa 

In the above examples, it becomes evident that the use of the 
development term produces an unacceptable amount of short-wave 


development. FNWF has a subroutine, SAG [i » which can be used 


to remove some or all perturbations. Itis of the form 
Ms 
A=At+tKOYA 


The problem remains as to what size to make k and how many 


times to use the smoother. 
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The smoother was tried with k = 0.125 and used three successive 
times on the D and H fields every six hours. The results of this heavy 


smoothing in addition to no sine function and K, at 0 is seenin Fig. 9. 


3 
Case IV 
Many of the systems have been smoothed out of existence and 
others have become unrepresentative of the verifying analysis. This 


example clearly shows that a low RMSE does not necessarily mean a 


good prognosis. 


mean deviation and range standard deviation 
Ve -7 lows 116.0(174, 92) Ban. | 
R 105 highs -29. 6(5, -77) 34,2 
Case V 


Maintaining the heavy smoothing, and setting K, at 1.0, starts to 


bring the prognosis pattern, Fig. 10, back to the analysis, but not 


enough. 
mean deviation and range standard deviation 
Pp 18. lows lltm 3( 142, 84) Coys 
JS ed l(Ohy) highs =2Onc(s, -a0) San 0 
Case VI 


Jiroie et akog, wba K, is 1.0 but is not diminished during the 24 hour 
period. As noted earlier, this would be too unstable to be computed, 
were it not for the use of a heavy smoother. 

Compared to Fig. 9, both highs and lows show much greater 
development, but this is attained with an excessive gain in short 


Waves. 
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mean deviation and range standard deviation 


P -8 lows 70. 9(117, 14) 47.2 
R 147 highs 8. 4(63, - 30) 300.5 
Case Vil 


Continuing the search for a compromise between development and 


smoothness, K, is again divided by 2 every 6 hours and the initial 


3 
value of K, is set at 2.0, which results in Fig. 12. While the central 


heights are closer to the analysis than Fig. 9 or 10, there is still too 


much short-wave development, although less than Fig. 11. 


mean deviation and range standard deviation 
P -8 lows 85. 5 (117, 55) 25.9 
RB. Gb2s highs 6, 2(53, -38) Bo.) 
Case VIII 


In Fig. 13, the smoothing coefficient is reduced to 1/16 vice 1/8 


while K., is 1.0. This is seen to produce development with acceptable 


Ss 
smoothness. 
mean deviation and range standard deviation 
P -8 lows 48. 9(69, 28) 14.7 
R 116 highs 4.0 (36, -21) teh ta 
Case IX 


Continuing in the same vein, the smoothing coefficient is reduced 
to 1/32 while K, remained 1.0, Fig. 14. Further development is noted 
but so is considerable short-wave formation. The compromise point 


has been passed. 
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mean deviation and range standard deviation 
P -8 lows . 3(26, - 22) ir Gs, 
R 140 highs 27. 6(65, 7) a6 
The above findings are now tested for a 48 hour period using a 
verification analysis of 1200Z 25 Jan. 1967, Fig. 15. The FNWF 
pillow and RMSE for this date were -16.0 and 173.2, respectively. 
The mean deviation, range and standard deviation of the lows were 
on 0(119,-230), 119 andiem the highs ~-Zeye eames gieend |]. 8. 
Case X 
Starting with the heavy smoother: k of 1/8, 3 scans every 6 hours, 
and K, at 0, results in Fig. 16. This is an extremely smooth chart 
with most baroclinic features missing. A misleading low RMSE is 
again encountered. 
mean deviation and range standard deviation 
1S 8 lows 117. 0(165, 81) OZ. 7 
me 143 highs -154(1 sample) = 
Case XI 


Continuing with the heavy smoother and changing K, to 1.0 produces 


3 
some slight development, Fig. 17, but still far short of the analysis. 
mean deviation and range standard deviation 
le 7 lows 101. 3(147, 66) 108.4 
R 142 highs -154(1 sample) - 
Case Xll 


Reducing the smoothing factor to 1/16 with K, of 0 improves the 


chart, Fig. 18, over the previous two results but still not enough. 
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mean deviation and range standard deviations 
1S 7 lows 67.0(112, 37) Os) 
R 148 highs -126(1 sample) = 
Case XIII 
With the smoothing coefficient at 1/16 and K, at 1.0, further 


development is observed in Fig. 19, but even more would be desirable. 


mean deviation and range standard deviation 
P ft lows 35, 4(84, - 5) Les, 
Ea 16] highs -67.0(-24, -110) 43 


An attempt was made at reducing k to 1/32 while leaving K, at 1.0. 
This variation produced a 30 hour prognosis before overflow occurred. 

Case XIV 

The final variation, Fig. 20, was the most satisfactory although 


not spectacular. Inthis variation, kis 1/16 but used twice instead of 


three times every six hours while K, is 1.0. 


mean deviation and range standard deviation 
P 8 lows - 55e4(- 145-23) eters: 
me 168 highs 5,0(62, -62) Bl. Ge 


The results of the above cases are summarized in Table III, 


3) | 
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7. Independent data results 

The last combination of factors, that is, K, of 1.0 reduced by one 
half every six hours and a smoothing coefficient of 1/16 applied to the 
height and thickness fields every six hours, was applied to independent 
data of seven consecutive 24-hour prognostic charts. 

These charts covered the period 0000Z 24 March to 0000Z 27 March. 
A comparison was made with full development term of 1.0, no develop- 
ment term and the FNWF product for these dates. 

Upon summarizing the results of the above tests, it was seen that 
a K, of 1.0 was probably too high and another group of seven charts 
was produced witha K, on-0.. 7, 

The results of all of the above are listed in Table IV. The pillow 
and RMSE for the variations on the applicable dates is presented in 
Tabie V. 

Table IV indicates a sharp reduction in the mean deviation of the 
lows from the analysis, however, the choice ofa K, of 0.7 was too 
much of a reduction from 1.0 in the case of the lows but was quite 
good for the highs. A K, of about 0. 8 would probably have produced 
better results for this period. 

Compared to FNWF, a K, of 0.7 results in a good reduction in the 


range of deviations for both highs and lows and a marked improvement 


for the standard deviation from the mean deviation for highs. 
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8. Conclusions 

The modification of the horizontal advection term by a function of 
latitude as used herein, generally improved the movement of troughs 
and ridges by reducing the lag between prognosis and analysis. 

It appears that a refinement of this function with a displacement 
of the maximum from about 20 degrees to about 45 degrees would be 
more rewarding. 

While compiling data for this study, there appeared to bea 
relation between the zonal index and the lag of the troughs and ridges. 

It is recommended that further study of this phenomena be conducted 
to arrive at a more exact relationship and how it could be applied to 
proper movement of troughs and ridges. 

It is readily seen that, as the name implies, the development term 
produces development in the form of deeper lows, higher highs and 
also, unfortunately, unwanted short waves. 

It results in distinct and selective changes in intensities of systems 
without materially changing the location of these systems. 

It showed marked improvement over the FNWF highs anda lower- 
ing of the range of errors in the lows. 

The inclusion of this term to some degree in an operational 
barotropic model would be beneficial. 

It is suggested that the most desirable approach would be to 
incorporate these empirical modifications in the FNWF 500-mb 
prediction equation while simultaneously carrying out the FNWF 1000-mb 


prediction hour by hour. The resulting H field could then be incorpor- 


Sif 


ated into the 500-mb forecasts each hour, which should be far more 
reliable than that produced by simple advection of thickness. 

A development term so produced could be made smaller initially, 
say K, = 0.25, but would not be reduced periodically since the thick- 
ness field would presumably stay in phase better with the 500-mb 
height field, 

This scheme should produce improved long range 500-mb and 


surface prognostic charts. 
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